Trans-Gulf migrants (TGMs) are species that migrate from Central and South America and land between Texas and Florida looking for areas to forage and breed. The specific time in which they land is imperative to their longevity, as arriving early could increase chances of finding a mate and increase time for mating, while arriving later would secure food availability. Population displacement is becoming a serious concern as migration patterns are being disrupted by the increasing severity of climate change, requiring the TGMs to alter the timing of their migration in order to access the same resources as they normally would. TGMs also rely on the environmental cues by stop over points throughout their trip to ensure that the migration will go smoothly, and the longer the migration trip, the less impact a shorter length shift in their departure date will have according to Carey (2009). To measure the effects of local weather conditions on TGM arrival time in Massachusetts, a program called eBird will be used for bird ecology data in tandem with meteorological data. In this project we will study the following five species: Magnolia warbler (Setophaga magnolia), Cape May warbler (Setophaga tigrina), Yellow-rumped warbler (Setophaga coronata), Palm warbler (Setophaga palmarum), and Black-throated Green warbler (Setophaga virens).
A Magnolia warbler.
Fig. 1.A. Records for Magnolia warbler in Massachussets, 2019.
A Cape May warbler.
Fig. 1.B. Records for Cape May warbler in Massachussets, 2019.
A Yellow-rumped warbler.
Fig. 1.C. Records for Yellow-rumped warbler in Massachussets, 2019.
A Palm warbler.
Fig. 1.D. Records for Palm warbler in Massachussets, 2019.
A Black-throated Green warbler.
Fig. 1.E. Records for Black-throated Green warbler in Massachussets, 2019.
After choosing what species of TGMs to study, occurence data from an rgbif function will be used to grab data for each species. NOAA’s NCDC API will be used to locate weather stations specific locations throughout the birds trip to Boston. Then, a function called meteo_pull_monitors() will be used to grab weather data from these stations. Finding the day where 25% of all the individuals of a species has arrived constitutes the “arrival time”, and the rest of the species that arrive will follow a logistic curve, modeled with a logistic model in R. After finding the weather data for these dates along the migration path of the species, a rolling mean function will help find the mean of the weather data over 2 weeks before each Julian day. A linear mixed-effect model will allow for analysis of many factors that may affect the relationship between arrival time and local weather. An ANOVA is run over the models to find a best fit model of the data.
## elevation mindate maxdate latitude name datacoverage station.type
## 1 65.5 1948-01-01 2022-12-11 30.68819 AL 1 GHCND
## 2 222.6 1939-01-01 2022-12-11 35.22254 NC 1 GHCND
## 3 3.2 1936-01-01 2022-12-11 42.36057 MA 1 GHCND
## id elevationUnit longitude longitude.1 latitude.1 migr.day
## 1 USW00013894 METERS -88.24598 1128508 -1510291.3 10
## 2 USW00013881 METERS -80.95433 1718307 -890812.8 5
## 3 USW00014739 METERS -71.00975 2320536 125815.9 0
Fig. 2. Weather sampling locations for this project.
## # A tibble: 6 × 46
## id date adpt aslp astp awbt awnd fmtm pgtm prcp rhav rhmn
## <chr> <date> <int> <int> <int> <int> <int> <int> <int> <dbl> <int> <int>
## 1 USW000… 2000-01-01 NA NA NA NA 35 1202 1341 0 NA NA
## 2 USW000… 2000-01-02 NA NA NA NA 45 1200 1157 0 NA NA
## 3 USW000… 2000-01-03 NA NA NA NA 68 1140 1418 0 NA NA
## 4 USW000… 2000-01-04 NA NA NA NA 50 40 1334 94 NA NA
## 5 USW000… 2000-01-05 NA NA NA NA 24 104 138 0 NA NA
## 6 USW000… 2000-01-06 NA NA NA NA 13 1006 1244 0 NA NA
## # … with 34 more variables: rhmx <int>, snow <dbl>, snwd <dbl>, tavg <dbl>,
## # tmax <dbl>, tmin <dbl>, tsun <int>, wdf2 <int>, wdf5 <int>, wesd <int>,
## # wsf2 <int>, wsf5 <int>, wt01 <int>, wt02 <int>, wt03 <int>, wt04 <int>,
## # wt05 <int>, wt06 <int>, wt07 <int>, wt08 <int>, wt09 <int>, wt11 <int>,
## # wt13 <int>, wt14 <int>, wt15 <int>, wt16 <int>, wt17 <int>, wt18 <int>,
## # wt21 <int>, wt22 <int>, wv03 <int>, psun <int>, wt19 <int>, wt10 <int>
First, a data set was created by collecting data for each of our species between 2000 and 2020 to allow to analyze eBird data. Then, another data set was created using NOAA’s NCDC API, collecting weather data in the between the years 2000 and 2020, in three different regions of the USA. The weather sampling locations are shown above, and were chosen to reflect the path of TGMs.
The plots above show that 25% of the population of Magnolia warblers arrived between the Julian days 130 and 136.
The plots above show that 25% of the population of Cape May warblers arrived between the Julian days 125 and 133.
The plots above show that 25% of the population of Yellow-rumped warblers arrived between the Julian days 114 and 121.
The plots above show that 25% of the population of Palm warblers arrived between the Julian days 103 and 112.
The plots above show that 25% of the population of Black-throated Green warblers arrived between the Julian days 125 and 131.
After analyzing both the eBird and weather data, linear mixed-effect modelling was conducted for the five species of birds.
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## tmin 0.2039 1 0.65160
## tmax 0.0049 1 0.94441
## wvec 0.4882 1 0.48473
## tmin:tmax 0.9260 1 0.33591
## tmin:wvec 1.0550 1 0.30437
## tmax:wvec 4.2667 1 0.03887 *
## tmin:tmax:wvec 0.6574 1 0.41749
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 0.7445 1 0.38824
## wk.tmax 0.4649 1 0.49532
## wk.wvec 1.2400 1 0.26547
## wk.tmin:wk.tmax 4.8273 1 0.02801 *
## wk.tmin:wk.wvec 1.0802 1 0.29866
## wk.tmax:wk.wvec 0.0987 1 0.75340
## wk.tmin:wk.tmax:wk.wvec 0.0632 1 0.80152
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
| (Intercept) | wk.tmax | wk.tmin | wk.wvec | wk.tmax:wk.tmin | wk.tmax:wk.wvec | wk.tmin:wk.wvec | wk.tmax:wk.tmin:wk.wvec | df | logLik | AICc | delta | weight | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 133.0531 | -0.0120013 | 0.0249803 | -0.0273876 | NA | NA | NA | NA | 6 | -125.2932 | 264.1713 | 0.00000 | 0.997795628 |
| 5 | 133.3764 | -0.0070455 | 0.0108047 | -0.1267917 | NA | NA | 0.0008826 | NA | 7 | -130.8749 | 277.9037 | 13.73244 | 0.001040110 |
| 3 | 133.1862 | -0.0081418 | 0.0145255 | -0.1434499 | NA | 0.000539 | NA | NA | 7 | -131.2900 | 278.7339 | 14.56266 | 0.000686751 |
| 2 | 124.1821 | 0.0182495 | 0.1322246 | -0.0247447 | -0.0003805 | NA | NA | NA | 7 | -131.6534 | 279.4607 | 15.28943 | 0.000477511 |
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 2.6329 1 0.1047
## wk.tmax 1.5240 1 0.2170
## wk.wvec 1.4580 1 0.2273
Based on the analysis shown above for the Magnolia warbler, a model including minimum and maximum temperature as well as wind vector explains the data best. The ANOVA test shows no significant predictors.
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## tmin 0.0134 1 0.9077
## tmax 0.0002 1 0.9880
## wvec 0.1326 1 0.7158
## tmin:tmax 4.7093 1 0.0300 *
## tmin:wvec 0.8355 1 0.3607
## tmax:wvec 0.0013 1 0.9714
## tmin:tmax:wvec 0.3551 1 0.5513
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 2.7879 1 0.09498 .
## wk.tmax 2.7641 1 0.09640 .
## wk.wvec 0.2455 1 0.62023
## wk.tmin:wk.tmax 0.5389 1 0.46288
## wk.tmin:wk.wvec 0.9111 1 0.33982
## wk.tmax:wk.wvec 0.1203 1 0.72871
## wk.tmin:wk.tmax:wk.wvec 0.0532 1 0.81755
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
| (Intercept) | wk.tmax | wk.tmin | wk.wvec | wk.tmax:wk.tmin | wk.tmax:wk.wvec | wk.tmin:wk.wvec | wk.tmax:wk.tmin:wk.wvec | df | logLik | AICc | delta | weight | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 129.7308 | -0.0222361 | 0.0370745 | -0.0198637 | NA | NA | NA | NA | 6 | -147.2659 | 308.1168 | 0.00000 | 0.9979509576 |
| 5 | 129.4959 | -0.0281771 | 0.0526071 | 0.0944574 | NA | NA | -0.0009794 | NA | 7 | -152.7535 | 321.6609 | 13.54417 | 0.0011429585 |
| 3 | 129.5709 | -0.0233099 | 0.0415797 | 0.0786820 | NA | -0.0004255 | NA | NA | 7 | -153.3383 | 322.8305 | 14.71372 | 0.0006368925 |
| 2 | 130.5758 | -0.0250536 | 0.0259247 | -0.0205305 | 3.94e-05 | NA | NA | NA | 7 | -154.1995 | 324.5528 | 16.43607 | 0.0002691914 |
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 2.6390 1 0.1043
## wk.tmax 2.1880 1 0.1391
## wk.wvec 0.2473 1 0.6190
Based on the analysis shown above for the Cape May warbler, a model including minimum and maximum temperature as well as wind vector explains the data best. The ANOVA test shows no significant predictors.
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## tmin 0.9037 1 0.3418
## tmax 0.6754 1 0.4112
## wvec 0.0733 1 0.7866
## tmin:tmax 0.0036 1 0.9522
## tmin:wvec 0.0283 1 0.8663
## tmax:wvec 0.0422 1 0.8372
## tmin:tmax:wvec 2.5256 1 0.1120
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 0.0224 1 0.8809
## wk.tmax 0.2524 1 0.6154
## wk.wvec 0.3964 1 0.5290
## wk.tmin:wk.tmax 0.0001 1 0.9906
## wk.tmin:wk.wvec 0.2760 1 0.5993
## wk.tmax:wk.wvec 0.3673 1 0.5445
## wk.tmin:wk.tmax:wk.wvec 0.0200 1 0.8874
| (Intercept) | wk.tmax | wk.tmin | wk.wvec | wk.tmax:wk.tmin | wk.tmax:wk.wvec | wk.tmin:wk.wvec | wk.tmax:wk.tmin:wk.wvec | df | logLik | AICc | delta | weight | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 118.0378 | -0.0058516 | 0.0044664 | -0.0198813 | NA | NA | NA | NA | 6 | -137.4599 | 288.5047 | 0.00000 | 0.9989817573 |
| 5 | 118.0398 | -0.0059926 | 0.0048242 | -0.0175320 | NA | NA | -2.44e-05 | NA | 7 | -143.8136 | 303.7810 | 15.27622 | 0.0004812467 |
| 3 | 118.0066 | -0.0051027 | 0.0026475 | -0.0505703 | NA | 0.0001442 | NA | NA | 7 | -144.0946 | 304.3431 | 15.83834 | 0.0003633341 |
| 2 | 118.2381 | -0.0065272 | 0.0009102 | -0.0199552 | 1.3e-05 | NA | NA | NA | 7 | -144.8328 | 305.8195 | 17.31476 | 0.0001736618 |
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 0.0791 1 0.7786
## wk.tmax 0.2923 1 0.5888
## wk.wvec 0.4214 1 0.5162
Based on the analysis shown above for the Yellow-rumped warbler, a model including minimum and maximum temperature as well as wind vector explains the data best. The ANOVA test shows no significant predictors.
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## tmin 0.0235 1 0.8781
## tmax 0.3636 1 0.5465
## wvec 0.6505 1 0.4199
## tmin:tmax 0.3784 1 0.5385
## tmin:wvec 0.2250 1 0.6353
## tmax:wvec 2.2826 1 0.1308
## tmin:tmax:wvec 0.9023 1 0.3422
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 3.6847 1 0.05491 .
## wk.tmax 10.6342 1 0.00111 **
## wk.wvec 0.2272 1 0.63364
## wk.tmin:wk.tmax 1.1488 1 0.28380
## wk.tmin:wk.wvec 0.1678 1 0.68204
## wk.tmax:wk.wvec 0.0028 1 0.95774
## wk.tmin:wk.tmax:wk.wvec 1.0208 1 0.31233
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
| (Intercept) | wk.tmax | wk.tmin | wk.wvec | wk.tmax:wk.tmin | wk.tmax:wk.wvec | wk.tmin:wk.wvec | wk.tmax:wk.tmin:wk.wvec | df | logLik | AICc | delta | weight | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 117.1578 | -0.0709318 | 0.0441983 | 0.0194340 | NA | NA | NA | NA | 6 | -138.4880 | 290.5609 | 0.00000 | 0.9987706585 |
| 5 | 116.9915 | -0.0687917 | 0.0400828 | -0.0037397 | NA | NA | 0.0003293 | NA | 7 | -144.6047 | 305.3632 | 14.80235 | 0.0006097850 |
| 3 | 117.1193 | -0.0701214 | 0.0418461 | -0.0195221 | NA | 0.0002159 | NA | NA | 7 | -144.9692 | 306.0923 | 15.53149 | 0.0004234920 |
| 2 | 115.2479 | -0.0632099 | 0.0841213 | 0.0186014 | -0.0001609 | NA | NA | NA | 7 | -145.7393 | 307.6325 | 17.07167 | 0.0001960645 |
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 4.8344 1 0.0278981 *
## wk.tmax 13.2222 1 0.0002766 ***
## wk.wvec 0.3294 1 0.5659997
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Based on the analysis shown above for the Palm warbler, a model including minimum and maximum temperature as well as wind vector explains the data best. The ANOVA test shows minimum and maximum temperature as significant predictors.
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## tmin 0.4563 1 0.4993
## tmax 1.1749 1 0.2784
## wvec 0.1405 1 0.7077
## tmin:tmax 0.3230 1 0.5698
## tmin:wvec 0.4894 1 0.4842
## tmax:wvec 0.0613 1 0.8044
## tmin:tmax:wvec 0.0759 1 0.7829
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 1.7645 1 0.1841
## wk.tmax 1.7308 1 0.1883
## wk.wvec 0.0955 1 0.7573
## wk.tmin:wk.tmax 0.5796 1 0.4465
## wk.tmin:wk.wvec 0.1317 1 0.7167
## wk.tmax:wk.wvec 0.3477 1 0.5554
## wk.tmin:wk.tmax:wk.wvec 0.8747 1 0.3497
| (Intercept) | wk.tmax | wk.tmin | wk.wvec | wk.tmax:wk.tmin | wk.tmax:wk.wvec | wk.tmin:wk.wvec | wk.tmax:wk.tmin:wk.wvec | df | logLik | AICc | delta | weight | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 128.5188 | -0.0115897 | 0.0172847 | 0.0034138 | NA | NA | NA | NA | 6 | -123.3140 | 260.2130 | 0.00000 | 0.9984621332 |
| 5 | 128.4319 | -0.0146353 | 0.0250088 | 0.0638277 | NA | NA | -0.0005165 | NA | 7 | -129.3160 | 274.7859 | 14.57293 | 0.0006836906 |
| 3 | 128.3984 | -0.0133190 | 0.0229514 | 0.1196528 | NA | -0.0005002 | NA | NA | 7 | -129.4846 | 275.1231 | 14.91012 | 0.0005776165 |
| 2 | 123.5982 | 0.0050947 | 0.0846200 | 0.0086033 | -0.000242 | NA | NA | NA | 7 | -130.2211 | 276.5961 | 16.38309 | 0.0002765598 |
## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: j.day
## Chisq Df Pr(>Chisq)
## wk.tmin 1.1974 1 0.2738
## wk.tmax 1.3751 1 0.2409
## wk.wvec 0.0161 1 0.8992
Based on the analysis shown above for the Black-throated Green warbler, a model including minimum and maximum temperature as well as wind vector explains the data best. The ANOVA test shows no significant predictors.
Weather can influence the arrival time for TGMs migrating to MA. In 2004, weather conditions were slightly above average for most of the year, including the Julian days between 120 and 150. That year is also said to be the fourth warmest year on record (US Department of Commerce). Our TGM data for three of our five species show a slightly different graphical display for those dates. Between the years 2000 and 2020, about 100% of our Magnolia Warblers, Cape May Warbler, and Black-Throated Green Warbler birds arrived in Massachusetts by Julian day 140. However, for all three species the 2004 data is an outlier. There appeared to be a delay for the three species in arriving in Massachusetts. There is a possibility that the 2004 higher than normal weather conditions globally confused these birds, delaying their migration patterns as weather does influence a bird’s arrival time. For the most part, minus the 2004 data, arrival time did not vary according to the temperature and wind variables along the migration route for the TGMs migrating to MA inspected. Over a span of 20 years, 5 TGMs: Setophaga magnolia, Setophaga tigrina, Setophaga coronata, Setophaga palmarum, Setophaga virens did not appear to have any significant changes within their migration patterns. Their arrival dates for 25% of their population varied only between a few days from year to year and their numbers always appeared to plateau around the same Julian day each year for the 20 recorded years.
Though climate change may play a role in the population status of TGMs arriving in MA during the spring, the data collected does not directly reflect climate change’s impact on arrival. All of our bird species data collected has their specific Julian day each year that is closest to 0.25 of the population arriving within a maximum of 8 days difference between the years 2000 and 2020. The year by Julian day plots for each bird species inspected do not show a clear pattern of the 25% arrival date being later in the year over time either. Instead, each TGM: Setophaga magnolia, Setophaga tigrina, Setophaga coronata, Setophaga palmarum, Setophaga virens plots look like a typical scatter plot with no visible potential line of regression. However, as previously mentioned, if years with weather similar to 2004 occur again, bird migration to the North may be delayed overall. While our data did not display a pattern of TGM migration being delayed the closer the year was to 2020, if temperatures are to rapidly increase due to global warming, bird migration to the North will slowly happen as less birds will arrive on their expected arrival dates as predicted by previous years.
In conclusion, our results do not convey a clear correlation between bird migration patterns and climate change. Instead, they support the idea that migration arrival dates vary slightly from the mean arrival date. A similar study mentioned how variation from the mean arrival date can be a source of more basic ecological changes that are not particular to climate change (MILLER-RUSHING et al. (2008)). Our specific findings pertaining to 2004 reflect that an increase in temperature directly affects bird migration. However, the data collected does not reflect an overall high enough North American temperature increase to display data which shows a pattern of bird migration arrival dates changing as according to climate change. Perhaps climate change is not at the detrimental point yet where global temperature increase directly modifies arrival times. Or perhaps, the data collected is not specific enough to showcase bird responses to climate change. A similar study looking at migratory birds and climate change suggests that perhaps, certain age groups or genders may respond differently to climate change (Knudsen et al. (2011)). It is possible that adults, having previously migrated to Massachusetts might arrive in Massachusetts on similar dates each year due to prior knowledge and experience. Younger birds with no prior migratory experience might have a slightly different variation to their arrival dates than their seniors which might reflect climate change’s influence on bird migration. While our results may not convey a clear correlation between bird migration patterns and climate change, they do indicate that high temperatures can delay bird migration. If global warming continues to drive climate change, our results suggest that bird migration to Massachusetts can be delayed as a whole.